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PRINCIPAL  NOMENCLATURE 


[Note:  Equation  numbers  are  given  below  where  it  may  help  to  identify  th 
point  of  first  introduction  of  a  symbol.] 


C  (  0 ) 

Curvature  of  body  surface 

c 

C(0) 

F(z) 

Complex  potential  corresponding  to  5(x,y) 

(24) 

fj  (z) 

Complex  potential  corresponding  to  <J>j  (x,y) 

(15) 

g 

Gravitational  constant 

gj  (z) 

Complex  potential,  a  part  of  f j (z)  (16) 

g(z) 

Complex- potential ,  a  part  of  F(z)  (25) 

H  (x) 

Free-surface  shape  from  exsct  solution  (5) 

H(X) 

Free-surface  shape  from  naive  expansion  (9) 

H  (x) 

Free-surface  shape  from  wavelike  solution 

$(x,y)  (21)  (31) 

hj  (z) 

Complex  potential,  a  part  of  f j  (z)  (16) 

h  (z) 

Complex  potential,  a  part  of  f(z)  (25) 

K 

d3R/d03  ,  evaluated  at  8=0  (37) 

k(z) 

[f q (z) ] (27) 

L 

Typical  body  dimension 

Pn(x) 

Function  defined  on  y  =  0  (13) 

R(Q^ 

Function  describing  body  shape  (1) 

Ro(0) 

Function  describing  shape  of  double  body  (34) 

r 

Radius  in  cylindrical  coordinates 

U 

Forward  speed  of  body  (or  speed  of  stream  past  fixed  body) 

x,y 

Cartesian  coordinates  (nondimensional) 

z 

x  +  iy 

Yn  (x»y) 

Re{gn (z) }  (46) 

e 

F2/gL  ,  the  small  parameter  of  the  problem 

Coefficient  of  en  in  the  naive  expansion  of  H(x)  (9) 

Phase  function,  -(l/e)/Zk(u)du  (28) 

z  1 

Angle  coordinate  in  cylindrical  coordinates 


nn(x) 

0(2) 

e 


Velocity  potential  for  the  exact  problem 
Velocity  potential  for  the  naive  expansion  (8) 
Velocity  potential  for  wavelike  motion  (20)  (30) 
Coefficient  of  en  in  naive  expansion  of  $\x,y) 
i?e{hn(x,y)  }  (47) 


I.  INTRODUCTION 

The  problem  of  ship-wave  generation  at  very  low  speeds  has  fundamental 
importance,  although  it  is  not  of  much  interest  to  naval  architects.  At  very 
low  speeds,  there  is  little  wave  generation  at  all,  and  so  low-speed  wave 
resistance  need  not  be  considered  in  the  design  of  a  ship.  However,  the  usual 
methods  of  predicting  wave  resistance  at  finite  speed  are  not  valid  as  speed 
approaches  zero,  that  is,  the  mathematical  solutions  on  which  these  methods 
are  based  are  not  uniformly  valid  at  low  speed.  So  there  is  a  question  how 
low  the  speed  must  be  for  the  standard  methods  to  break  down. 

There  is  another  practical  reason  for  investigating  the  low-speed  problem: 
Inui  and  Kajitani  (1977)  found  that  the  propagation  of  ship-generated  waves 
can  be  predicted  more  accurately  if  the  waves  are  considered  to  propagate  on 
the  nonuniform  stream  around  a  double  body  (the  actual  part  of  the  body  below 
the  undisturbed  free  surface  plus  its  mirror  image)  rather  than  on  a  uniform 
stream.  A  low-speed  analysis  leads  naturally  to  such  an  approach.  It  may  be 
that  "low  speed"  to  a  mathematician  may  not  be  very  low  at  all  to  a  naval 
architect. 

Over  a  period  of  decades,  numerous  authors  have  discussed  the  low-speed 
limit  of  solutions  obtained  by  standard  perturbation  techniques.  Such  solutions 
are  based  on  the  assumption  that  all  perturbation  velocities  are  asymptotically 
small  in  magnitude  compared  to  the  forward  speed  U  of  the  body  (or,  if  the 
reference  frame  is  fixed  to  the  body,  compared  to  the  speed  U  of  the  streaming 
flow  past  the  body) .  The  implications  of  this  assumption  become  unclear  if  then 
one  lets  U  -*■  0  . 

Salvesen  (1969)  showed  for  a  particular  case  that  this  is  a  singular  limit. 
Initially,  he  considered  the  forward  speed  as  fixed.  Then  he  expanded  the  prob¬ 
lem  statement  in  terms  of  a  small-disturbance  parameter,  worked  out  the  explicit 
solution  to  three  terms,  and  finally  examined  the  behavior  of  the  expansion  as 
U  ■*  0  .  At  moderate  speed,  the  second  and  third  terms  provided  small  correc¬ 
tions  to  the  linear  solution.  Then,  as  U  decreased,  the  second-order  term 


2 


became  larger  and  larger  until  it  dominated  the  first-order  term.  At  still 
lower  speed,  the  third-order  term  began  to  dominate  the  .-econd-order  term. 

It  appeared  that  the  problem  was  becoming  more  and  more  nonlinear  as  speed 
decreased  more  and  more.  However,  from  a  physical  point  of  view,  it  is  clear 
that  the  fluid  disturbance  becomes  smaller  and  smaller  as  speed  decreases. 

This  apparent  paradox  is  a  result  f  nonuniformity  of  the  expansion  as  U  ■*  0  . 

Ogilvie  (1968)  proposed  a  method  of  solution  that  was  expected  to  remove 
the  nonuniformity.  He  reasoned  as  follows:  At  very  low  speed,  the  free  sur¬ 
face  is  hardly  disturbed  at  all,  and  so  a  first  approximation  should  be  equiva¬ 
lent  to  the  double-body  flow.  In  such  a  flow,  if  the  fluid  is  ideal,  the 
magnitude  of  the  perturbation  velocity  is  everywhere  proportional  to  U  .  So 
one  should  not  assume  a  priori  that  the  perturbation  velocities  vanish  much 
more  rapidly  than  U  in  the  limit.  Rather,  one  should  take  the  double-body 
flow  as  the  first  approximation.  Then  perhaps  it  is  reasonable  to  assume  that 
free-surface  effects  are  small  perturbations  of  this  nonuniform  base  flow. 

This  was  a  new  assumption,  which  Ogilvie  implemented  on  a  highly  pragmatic 
basis :  He  argued  that  there  was  little  chance  of  solving  the  wave  problem 
unless  it  were  a  linear  problem,  and  so  he  tried  to  determine  under  what  condi¬ 
tions  this  would  in  fact  be  the  case,  and  he  obtained  a  partial  solution  for 
the  problem  of  a  fully  submerged  two-dimensional  hydrofoil  (Salvesen's  problem) 

Effectively,  Ogilvie  divided  the  flow  field  into  two  parts,  the  two  being 
characterized  by  vastly  different  length  scales.  In  the  double-body  flow,  the 
fluid  motion  can  be  described  on  a  scale  comparable  to  body  dimensions.  The 
wave  motion,  however,  is  described  appropriately  on  a  scale  comparable  to  wave¬ 
length.  A  suggestion  of  what  this  scale  ought  to  be  can  be  obtained  from 
classical  (finite-speed)  wave  theory:  Waves  on  a  stream  of  speed  U  typically 
have  wavenumber  <  =  g/U2  and  the  corresponding  wavelength  is  X  =  2tt/k  = 
2nU2/g  ,  where  g  is  the  gravitational  constant.  It  seems  appropriate  then  to 
introduce  a  small  parameter  c  =  F2  =  U2/gL  ,  where  L  is  a  typical  body  dimen 
sion  and  F  is,  of  course,  a  Froude  number.  In  the  wave  part  of  the  solution, 
differentiation  of  a  field  variable  has  an  order-of-magnitude  effect  equivalent 
to  multiplication  by  the  wav'numer,  which  is  0(e  *)  .  This  fact  was  used 
systematically  to  order  terms  in  the  several  conditions  that  constitute  the 
statement  of  the  wave  problem. 


This  division  of  the  problem  into  two  distinct  [■arts  is  based  on  a 
plausible  argument,  but  several  fundamental  questions  remain  unanswered: 

(1)  One  can  develop  a  perturbation  expansion  of  the  solution  without 
introducing  the  concept  of  the  rapidly  varying  wave  field.  The  expansion  thus 
developed  has  been  called  the  "naive  expansion."  It  appears  on  first  sight  to 
satisfy  the  formal  conditions  of  the  problem.  However,  it  cannot  possibly 
represent  a  wave  motion  (with  short  waves),  and  so  it  is  incomplete  or  possi¬ 
bly  quite  incorrect.  How  can  one  show  formally  that  it  is  incomplete  if  it 
satisfies  the  partial  differential  equation  and  all  of  the  boundary  condition? 
How  does  one  determine  the  relationship  between  it  and  a  truly  wave  like 
solution? 

(2)  Dagan  (1972)  follows  much  the  same  reasoning  as  Ogilvie,  but  he 
argues  (persuasively)  that  the  base  flow  should  correspond  to  the  two-term 
naive  expansion. 

(3)  Keller  (1979)  developed  a  ray  theory  based  on  the  same  fundamental 
concept.  Like  Dagan,  he  used  the  two-term  naive  expansion  for  the  base  flow. 
However,  he  argued  that  Ogilvie 's  linear  free-surface  condition  on  the  wave 
solution  is  defective:  It  contains  terms  on  one  side  that  are  rapidly  vary¬ 
ing,  but  the  terms  on  the  other  side  of  the  equation  are  not  rapidly  varying, 
and  so,  as  Keller  pointed  out,  they  cannot  be  equal  (unless  both  sides  equal 
zero,  which  is  not  an  acceptable  resolution  of  the  difficulty) .  So  Keller 
set  up  the  wave  problem  in  such  a  way  that  the  free-surface  condition  would 
be  homogeneous  in  the  rapidly  varying  field  quantities.  This  then  leads  to 
difficulties  in  predicting  how  waves  are  generated.  The  only  mechanism  (from 
a  mathematical  point  of  view)  must  exist  at  the  body  surface.  However,  a  ray 
theory  cannot  be  used  to  predict  what  happens  near  a  body  in  the  presence  of 
the  free  surface,  since  the  governing  partial  differential  equation  is  ellip¬ 
tic,  and  the  pseudo-hyperbolic  character  of  the  ray-theory  solution  is  not 
valid  there.  The  ray  theory  can  be  used  only  to  predict  how  water  waves  will 
propagate  once  they  have  been  generated. 

(4)  It  is  not  clear  whether  Ogilvie 's  approach  resolved  the  apparent 
paradox  in  Salvesen's  perturbation  expansion. 


In  this  paper,  we  shall  provide  some  answe: /.  to  the.  -  nie  d;ons,  although 
the  actual  solution  is  presented  is  a  subsequent  paper,  Chen  and  Ogilvie  (  1j8 J ) . 
In  particular,  we  shall  shew  the  following: 

(1)  The  naive  expansion  in  fact  does  not  satisfy  all  cf  the  boundary  con¬ 
ditions  of  the  problem. 

(2)  In  showing  that  the  boundary  conditions  are  not  completely  satisfied 
by  the  naive  expansion,  we  obtain  a  nonhomogeneous  boundary  condition  that  will 
have  to  be  satisfied  by  the  wave  like  part  of  the  solution. 

(3)  Asymptotically  as  0  ,  Ogilvie's  solution  does  not  represent  a 

wavelike  motion.  We  now  believe  that,  in  a  strict  asymptotic  sense,  Uie  wave 
resistance  of  a  submerged  body  is  small  of  exponential  order  in  this  limit. 

We  have  not  been  able  to  obtain  a  relationship  between  the  expansion  for  e  ->■  0 
and  Salvesen's  expansion. 

(4)  The  double-body  flow  is  not  analytic  at  the  intersection  of  the  body 
and  the  undisturbed  free  surface.  The  singularity  is  very  weak,  but  subsequent 
terms  in  the  naive  expansion  exhibit  singularities  there  that  are  less  weak. 

Following  Keller  (1979) ,  we  conclude  that  the  free-surface  condition  on 
the  wave  part  of  the  solution  must  indeed  be  homogeneous.  Wave  generation  is 
accomplished  through  the  body  boundary  condition.  Baba  and  Takekuma  (1975) 
and  Maruo  and  Fukazawa  (1979)  followed  Ogilvie's  approach,  in  which  the  waves 
arise  from  a  nonhomogeneous  free-surface  condition;  their  results  appear  to 
be  questionable. 

Two  more  points  should  be  noted:  (i)  Very  short  waves  are  strongly  affected 
by  surface  tension,  which  we  neglect.  But  wa  are  interested  ultimately  in  the 
problem  of  ship-generated  waves,  and  even  very  "short"  waves  in  that  case  should 
be  long  enough  not  to  be  subject  to  serious  alterations  by  surface  tension.  (ii 
In  the  submerged-body  problem,  our  conclusions  may  not  be  val’d  if  the  submer¬ 
gence  is  very  small  (comparable  to  wavelength) . 


II.  PROBLEM  FORMULATION 

We  consider  a  slow  streaming  flow  past  a  two-dimensional  body.  Let  the 
body  surface  be  given  by  an  equation  of  the  form 

r  =  R (0 )  ,  (1) 

where  r  =  (x“+y2)  1//2  and  0  is  the  polar  angle  measured  counterclockwise 
from  the  x  axis.  The  stream  has  speed  U  in  the  positive  x  direction. 

The  undisturbed  free  surface  coincides  with  the  x  axis. 

All  lengths  are  nondimensionalized  by  a  length  L  ,  which  may  be  taken  as 
any  characteristic  body  dimension.  We  seek  an  asymptotic  solution  expressed 
in  terms  of  the  parameter 

£  =  F2  =  U2/gL  ,  (2) 

where  F  is  the  Froude  number  and  g  is  theacceleration  of  gravitv.  The 
velocity  potential  is  taken  as  LU4(x,y)  ,  and  so  4>  is  also  nondimensional . 
The  shape  of  the  free  surface  will  be  given  as  y  =  cH(x)  . 

The  potential  satisfies  the  Laplace  equation, 

[L]  $xx  +  <*>yy  =  0  in  the  fluid  domain,  (3) 

the  body  boundary  condition, 
r  1  34> 

[B]  —  =0  on  the  body,  (4) 

dn 

the  dynamic  free-surface  condition, 

[H]  H  (x)  =  J  *y}|  ,  (5) 

'y=eH(x) 

the  free-surface  condition  combining  [H]  with  the  kinematic  condition: 

[F]  4>y  +  +  2l>x1>y4’Xy  +  4>y<l>yy}  =  0  on  y  =  cH(x)  ,  (6) 

and  a  radiation  condition, 

[R]  |  $  -  x|  -*  0  as  x  -e»  and/or  y  -<*>  .  (7) 


The  naive  expansion  exhibits  one  kind  of  approach  to  .--oiving  such  a  pr  >bj.. 
Assume  that  ^>(x,y)  and  H(x)  can  be  expanded  i.i  series  in  t.  ; 


N 

<Mx,y) 

=  J>(x,y;e) 

n  =  c 

N 

(d) 

H  (x) 

=  II(x;e) 

'  1  tnnn(x)  . 

n=  0 

(9) 

When  these  expansions  are  simply  substituted  into  the  problem  statement,  one 
obtains  the  following  set  of  problems: 


[B] 

—  —  =0  on  the  body : 

jn  J 

(11) 

[H] 

Vx)  =  I  i1"  ^Oxl  ' 

(12a) 

Vx)  =  -  4>0x4l}x  , 

(12b) 

n2(x)  = 

“  *0X*2X  “  I  "  I  ”  ^0*1  t<*i0x+lxy  +  ♦ly*0yy2 

~  o'  I1  "  ‘J’OjJ  ^Ox^Oxyy  +  • 

(12c) 

all  to  be  satisfied  on  y  =  0  ; 


The  first  term  in  the  expansion,  <j>0(x,y)  ,  is  the  solution  of  the  rigid- 

wall  problem,  in  which  the  free  surface  is  replaced  by  a  wall.  We  sometimes 
refer  to  this  problem  also  as  the  "double-body"  problem.  Notwithstanding  the 
interpretation  in  terms  of  a  hypothetical  rigid  wall  at  y  =  0  ,  there  is  a 
corresponding  free-surface  disturbance,  given  by  (12a).  However,  it  is  not  a 


wavelike  disturbance.  In  fact,  no  matter  how  far  this  method  of  solution  is 
pursued,  there  will  never  be  any  real  waves. 


Clearly,  this  is  not  the  correct  solution,  since  we  expect  to  find  waves 
behind  the  body.  Nevertheless,  the  expansion  (8)  will  have  an  important  role 
later,  and  so  we  note  that,  in  principle,  it  can  be  solved  completely.  The 
4>q  problem  can  be  solved  by  standard  numerical  methods.  To  proceed  further, 
let 


<t> j  (x,y)  =  Re{ fj(z)}  ,  (15) 

where  f j (z)  is  a  function  of  the  complex  variable  z  =  x  + iy  .  Also,  set 

f  j  (z)  =  gj  (z)  +  hj  (z)  ,  j  >  0  ,  (16) 


where 

ds  pj (s) 

s  -  z  ' 

F 

and  F  is  the  part  of  the  x  axis  outside  of  the  body.  Then  we  have 


h-:  (z)  =  - 

J  IT 


(17) 


hj(z)  -  7 


f  ds  pj (s) 
s  -  z 


-  ip 


j(x)  +7- 


i  [  ds  Pj(s) 


S  -  X 


(18) 


f  yf° 

the  last  integral  being  a  principal-value  integral.  When  these  results  are 
substituted  into  (13b)  and  (13c) ,  we  find  that 


Im{g\ (x-iO) }  =  0  .  (19) 

Thus  gj  (z)  satisfies  conditions  just  like  fQ(z)  »  the  body  boundary  condi¬ 
tion  being  written  so  that  the  combination  9j+hj  satisfies  (11) .  The  numerical 
method  used  to  solve  for  fQ(z)  can  ^  use(l  f°r  all  of  the  gj  (z)  . 

An  alternative  method  of  obtaining  hj (z)  will  be  presented  later. 

The  expansions  (8)  and  (9)  are  not  complete.  It  will  be  necessary  to 
augment  them  with  terms  that  represent  real  waves  and  to  obtain  appropriate 
conditions  to  be  satisfied.  An  attempt  to  do  this  was  presented  by  Ogilvie 
(1968) ,  and  several  investigators  have  followed  his  approach.  Unfortunately, 
this  does  not  lead  to  the  desired  results.  It  is  worthwhile  to  consider  why, 
which  we  do  in  the  next  section. 


III.  THE  ANALYSIS  OF  OGILVIE  (1968) 


II 


« 


fi 
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Ogilvie  assumed  that  the  solution  consisted  of  (i)  the  rigid-wall  solution 
and  (ii)  a  wave  term: 

lMx,y)  =  4>0  (x,y)  +  e2<f>  (x,y;  e)  ,  (22) 

H  (x)  =  n0(x)  +  cH(x;e)  .  (21) 

The  new  terms  were  assumed  to  have  the  property  that 

$x  ,  5y  =  0($/e)  ;  H'  =  0(H/e)  .  (22) 

This  assumption  corresponds  to  the  expectation  that  the  waves  will  be  very 
short,  with  wavelength  that  is  0(e)  (and  thus  wavenumber  that  is  0(e_1)  ). 

When  (20)  and  (21)  are  substituted  into  the  free-surface  conditions,  the 
free-surface  condition  is  found  to  be: 

e2$y+ e34»02x$xx  =  e  •  (23) 

In  this  condition,  <}> Q  is  to  be  evaluated  on  y  =  0  and  $  on  y  =  enQ  (x)  . 
In  principle,  it  is  wrong  to  evaluate  $  on  y  =  0  ,  but  Ogilvie  (and  others 
following  him)  have  devised  arguments  indicating  that  perhaps  this  can  be  done 
anyway . 


For  the  case  of  a  submerged  body  (the  only  case  considered  by  Ogilvie) , 
one  can  construct  a  partial  solution,  as  shown  by  Ogilvie.  Let 


^(x,y;e)  =  Re{ F(z)}  , 


(24) 


F(z)  =  g(z)  +  h  (z)  , 


(25) 


with  h(z)  given  by: 


h'(z)  = 


Trie 


,  Z 


ds  p J (s) 


dc  kU)  itG(z)-eU)] 

.  e 

s  -  e 


(26) 


and 


k ( z )  =  [fg(z)]'2  , 


(27) 


1 

"  e  .. 


k(u)  du 


61 


-  ft  - 


0  (z) 


(28) 


The  value  of  z ^  is  arbitrary,  since  it  vanishes  when  the  difference, 

[0(z)  -0(c)]  is  introduced.  When  h' (z)  is  used  in  (23),  it  is  found  that  it 
satisfies  the  nonhomogeneous  condition,  leaving  only  a  homogeneous  free-surface 
condition  to  be  satisfied  by  g(z)  .  Of  course,  h(z)  does  not  satisfy  the 
body  boundary  condition  in  general,  and  g(z)  must  be  determined  so  that  the 
complete  solution  does  satisfy  that  condition. 

The  partial  solution  h(z)  was  constructed  in  much  the  same  way  that  one 
finds  solutions  of  problems  with  pressure  distributions  imposed  on  the  free 
surface  (see  Wehausen  and  Laitone  (1960) ) .  So  one  might  assume  that  it  repre¬ 
sents  a  true  wave  motion  that  decays  rapidly  with  depth.  If  that  were  true, 
the  term  g(z)  would  not  be  really  necessary,  since  h(z)  represents  extremely 
short  waves  (if  it  represents  waves  at  all).  In  such  a  case,  we  would  have 
found  the  solution. 


In  fact,  if  one  were  to  compute  the  solution  represented  by  h(z)  ,  it 
would  undoubtedly  exhibit  a  wavelike  behavior.  But  this  is  not  true  in  an  asymp¬ 
totic  sense,  as  we  now  show.  First,  we  note  that  h' (z)  appears  superficially 
to  be  0(e-1)  (see  (26)).  Now  integrate  by  parts,  first  with  respect  to  c  and 
then  with  respect  to  s  .  The  result  is  the  following: 


h’  (z) 


1  f  ds  pj(s) 


iQ(z) 


TT  I  S  -  Z 
i.  00 


-i©U) 


ds  pj (s) 

s  —  C 


The  first  term  is  clearly  0(1)  and  it  does  not  represent  a  wave  motion.  The 
second  term  appears  to  be  of  the  same  order  of  magnitude,  but  another  integration 
by  parts  (with  respect  to  both  c  and  s  )  leads  to  the  result: 


h'  (z) 


I 
*  . 


1  f°ds  pj  (s) 


,  e 

1£ 

k(z) 


ds  P'i(s)  iee 


iO(z) 


-i0(C) 


to 

ds 

’J_s  ■ 


ds  fp'J'ls)  p'i  (s)k'(O 


Cl  k(t) 


k2U) 


So  the  second  term,  also  not  wavelike,  is  0(e)  .  The  process  can  clearly  be 
continued  indefinitely,  and  we  see  that  we  obtain  a  series  in  increasing  powers 
of  e  ,  none  of  the  terms  representing  a  wave  motion. 

The  solution  in  (26)  corresponds  closely  to  the  naive  expansion.  The  first 
term  in  (29)  is  identical  to  h’(z)  in  (18).  The  following  terms  in  (29)  could 
be  obtained  from  a  naive  expansion  of  the  linear  problem  stated  in  (23) . 


-  In  - 


Thus  Ogilvie's  (1968)  approach  does  not  solve  t.ne  p roule.n  of  determining 
the  wave  motion  at  very  low  speed.  The  results  ol  evaluating  (10)  numerically 
are  not  known,  and  they  just  might  give  reasonable  predictions  in  fne  actual 
physical  situation,  although  they  would  be  inconsistent  from  a  perturbation- 
theory  point  of  view.  It  also  appears  that  reasonable  results  might  be 
obtained  if  the  body  were  sur face-piercing,  for  then  the  integrations  by  ports 
would  not  yield  the  trivial  result  found  in  (29).  But  that  is  not  really  :..;e 
direction  in  which  to  look  for  a  valid  abroach. 

A  fundamental  error  in  the  above  analysis  appears  to  have  been  made  in  the 


assumption  in  (20)  and 

(21)  . 

une  should  write: 

<t  (x, y) 

=  $(x,y?t)  +  f.2Hx,y;£)  , 

(30) 

H  (x) 

=  H (x;  e)  +  eH(x;c)  . 

(31) 

That  is,  the  wave  term 

must  be 

added  to  the  complete  naive  expansion 

(or  to 

at  least  two  terms  in  it).  When  (30)  and  (31)  are  substituted  into  the  free- 
surface  conditions,  one  obtains: 

e2$y  +  e3$^xx  =  0  .  (32) 

That  is,  the  free-surface  condition  is  homogeneous.  The  right-hand  side  of  (23) 

lias  been  absorbed  into  the  problem  for  4> ^  in  the  naive  expansion  (see  (13b)). 

Keller  (1979)  argued  this  for  different  reasons.  It  may  also  be  noted  that 

-2 

Dagan  (1972)  insisted  that  the  coefficient  in  (32)  (or  (23))  should  be  $x 

2 

and  not  vq^  ,  again  for  different  reasons.  We  agree  with  Keller  and  with 
Dagan. 

The  implication  of  (32)  is  that  there  is  no  "equivalent  pressure  distribu¬ 
tion"  on  the  free  surface  causing  the  ere* '.ion  of  waves.  In  the  case  of  the  sub¬ 
merged  body,  there  appears  to  be  no  mechanism  that  can  lead  to  the  creation  of 
waves  —  a  conclusion  that  is  compatible  with  the  result  expressed  in  (29) .  For 
the  case  of  a  surface-piercing  body,  we  may  expect  to  find  a  nonhomogeneous  con¬ 
dition  in  the  neighbor hood  of  the  intersection  of  the  body  and  the  free  surface. 
This  we  shall  show  explicitly  in  the  next  section. 


IV.  A  NONHOMOGENEOUS  CONDITION  FOR  THE  WAVE  POTENTIAL 

As  already  noted,  it  appears  that  the  "naive"  expansion,  (8)  and  (9) ,  satis¬ 
fies  the  boundary  conditions  of  the  problem,  leaving  a  homogeneous  problem  for 
the  wave  potential.  This  is  an  unsatisfactory  situation,  since  it  is  not  clear 
then  how  the  waves  ever  come  into  existence.  We  need  a  nonhomogeneous  boundary 
condition  to  impose  on  the  wave  potential. 

A  closer  examination  of  the  naive  expansion  will  show  that  it  does  not  com¬ 
pletely  satisfy  the  body  boundary  condition  in  the  case  of  a  surface-piercing 
body,  it  is  the  purpose  of  this  section  to  show  this  and  to  derive  the  nonhomo¬ 
geneity  in  the  wave  problem. 

The  clue  to  accomplishing  this  task  can  be  found  in  (13a)  :  If  we  find  a 
function  <j>Q  that  satisfies  the  Laplace  equation  in  the  lower  half-space  out¬ 
side  of  the  body,  that  function  can  be  continued  analytically  into  the  upper 
half-space  because  of  its  property  in  (13a) .  We  suppose  that  this  function 

satisfies  the  body  boundary  condition  in  y  <  0  .  There  is  no  guarantee  that 

this  function  also  satisfies  the  body  boundary  condition  in  y  >  0  .  And  it 

must  be  noted  that  part  of  the  wetted  surface  of  the  body  does  in  fact  lie  in 

the  upper  half-space,  since,  from  (5),  H(x)  =  1/2  (and  thus  y  =  e/2  )  at  a 
stagnation  point.  Similar  conclusions  can  be  drawn  from  (13b),  (13c) ,  etc., 
although  the  situation  is  not  so  obvious. 

A.  THE  $0  PROBLEM 

From  (13a),  we  can  write: 

<J>Q(x,y)  =  <j>0(x,-y)  .  (33) 

We  assume  that  <j>Q  satisfies  the  Laplace  equation  in  the  lower  half-space  (out¬ 
side  of  the  body  and  that  it  satisfies  the  body  boundary  condition  there,  (11). 
Our  purpose  no-  s  to  compute  34> Q/3n  on  that  part  of  the  body  for  which  y  >  0  . 

We  assume  for  now  that  the  function  defining  the  body,  R(0)  (see  (1)), 

has  at  least  three  continuous  derivatives  in  a  range  | 9 |  <  e/2  .  The  most 


11  - 


12 


interesting  case  is  that  d3R/d63  ?  0  at  6=0.  We  assume  that  the  body  sur¬ 
face  is  vertical  at  6=0,  which  implies  that  dR/dG  =  0  at  0  =  J  . 


As  already  noted,  represents  the  streaming  flow  around  the  douole  body. 

In  -tt  <  G  <  0  ,  (1)  describes  both  the  actual  body  and  the  double  boa.  .  We 

introduce  the  following  description  of  the  double  body,  applicable  in  bo.:h  6  <  0 
and  6  >  0  : 


r  =  R0(e)  = 


R(0)  -r  <  3  v  0  , 

R(-0)  ,  0  <  6  <  it  . 


(34) 


In  a  neighborhood  of  0  =  0  ,  we  can  expand  R(G)  and  RQ(6) 

R(0)  =  R(0)  +  02  R"  (0)  +  — -  63  K  +  ...  , 

^  C.2  Dn  / r\\  _  !  n  I  3 


Rq  ( 0)  =  R  (0)  +  —  9^  Ru  (0)  -  —  !  0  I  K  +  •  •  •  , 


where 


K  = 


d  3  R 

de3 


(35) 

(36) 

(37) 


6=0 


We  now  compute  34>Q/3n  on  r  =  R(6)  .  This  quantity  vanishes  for  6  <  0  ,  but 
not  otherwise.  For  the  moment,  it  is  convenient  to  write  <j>0  =  4>Q  (r,  0)  .  In 
polar  coordinates,  we  can  write  the  normal  derivative:* 


(38) 


3<j)0  _  -R(6)  <j)0r  +  [R*  (8)/R(6)  ] 

3n  {  [R(0))2+  [R’  ( 6)  1  2 } 1/2 


We  now  evaluate  (j>Qr  (R(0)  ,  0)  in  terms  of  values  on  (RQ(0),0)  .  For  example, 
4>Or(R(0),0)  =  *0  (RO(0),0)  +  [  R  ( 8 )  -  Rq  (8)  ]  i)>orr  (Ro  +  •••  •  (39) 

From  (35)  and  (36) ,  we  have 

1 

0  >  0  . 


R(  0)  -  Ro(0)  =  ~  K  6 3  + 


Carrying  out  details  in  this  way,  we  find  that 


Hjl 

dn 


r=R(0,  tRo(0)i2  r  TUe 


024)Ofl(Ro(0)  ,0)  +  70'  <f>0 


00 


+  O(04) 


(40) 


(41) 


We  note  also  that 


We  take  the  unit  normal  outward  from  the  fluid. 


Rq(0)  =  R0(O)  +  0(0*)  , 

<t>Oe(Ro(0)  ,0)  «  <t>O0(RO(O)fO)  +  e{Ro(O)4>or0  +  4>oe0}  +  0(63) 

=  0*O0e(RO(O)  +  °(03)  * 

the  last  result  following  from  the  facts  that  ^Og (R0 (°) »0)  =  0  and  Rg (0)  =  0 
Thus  (41)  simplifies: 


1  r=R(0) 


VTTitrp  e3»oee(Ro(0)'0)  *  0161,1  "  ole3) 


This  quantity  will  have  to  be  canceled  by  a  contribution  from  the  wave  potential. 

If  it  happens  that  K  =  0  ,  the  leading-order  term  in  (42)  will  be  0(e5)  , 
provided  that  d5R/d05  /  0  at  9=0. 

In  a  formal  sense,  one  can  claim  that  4>q  does  indeed  satisfy  the  body 
boundary  condition,  even  for  y  >  0  ,  since  <J>Q  is  simply  the  first  term  in  an 
asymptotic  solution,  and  the  body  boundary  condition  requires  only  that  3<(>0/9n 
=  o(l)  on  r  =  R( 0)  .  This  is  a  legitimate  point  of  view,  of  course.  If  we 
follow  it,  we  would  then  expect  the  quantity  in  (42)  to  show  up  in  the  statement 
of  a  higher-order  problem,  and  we  might  suppose  that  we  could  then  cancel  it  with 
an  appropriate  contribution  from  the  higher-order  solution.  However,  it  does  not 
work  out  this  way.  At  each  stage  of  the  solution  of  this  problem,  we  shall  find 
another  quantity  more  or  less  like  that  in  (42),  always  in  violation  of  the  body 
boundary  condition  for  0  <  y  <  e/2  . 

Another  special  case  is  possibly  of  comparable  interest:  Let  the  shape  in 
y  <  0  be  described  as  above,  but  assume  that  the  body  is  continued  above  y  =  0 
by  vertical  walls.  Then  we  can  write  for  the  body  geometry: 


r  =  R(  0 )  = 


Ro(0) 


0  <  0  , 


R0  (0)  sec  0  ,  9  >  0 


In  this  case,  the  body  curvature  is  not  generally  continuous  at  9 
the  general  formula  for  curvature, 


C  ( 0)  = 


[R(9)]2  -  R  (  9 )  R1'  ( 9 )  -  2  [R1  (9)  1  2 


{ [R( 9) ] 2  +  [ R" ( 9 ) ] 2  } 


2  \  3/2 


we  have  for  9  t  0  : 


C(-0) 


R(0)  -  R" (-0) 
R2  (0) 


(44) 


whereas  for  910  we  have: 


C(+0)  =  0  . 

If  these  happen  to  be  equal,  the  preceding  analysis  covers  the  case.  So  lot 
us  now  assume  that  the  limiting  curvature  from  below  is  not  zero.  The  analysis 
is  still  similar  to  the  preceding,  and  we  find  that,  for  9  >  0  , 

=  C (-0)  e2  |j[R(0)]  ?4>0rr(R0(O5,O)  -  *0  }  +  O(03)  .  (45) 

r=R ( 0)  rr 

So  the  error  in  satisfying  the  body  boundary  condition  above  y=0  is  0(e2) 

in  this  case. 

In  deriving  some  of  the  above  results,  we  implied  that  4>o(r,0)  can  be 
expanded  in  Taylor  series.  This  is  not  really  true  at  the  stagnation  points, 
where  in  fact  this  potential  is  not  analytic  unless  the  function  RQ(0)  is 
also  analytic.  However,  the  singularity  is  extremely  weak,  and  the  expansions 
used  above  are  valid  as  far  as  they  have  been  carried.  This  point  is  discussed 
further  in  the  Appendix. 


lio 
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B.  THE  4>  PROBLEM 

In  the  4^  problem,  there  is  an  error  in  satisfying  the  body  boundary  con' 
dition  just  like  the  error  described  above  in  the  4>q  problem.  Since  we  must 
multiply  <J>i  by  e  in  the  expansion  of  4>(x,y;e)  ,  this  error  is  0(e4)  (for 
the  case  of  continuous  curvature  at  0=0). 

There  is  another  error  in  the  4>j  problem.  Because  of  the  nonhomogeneous 
free-surface  condition  on  4>j  <  (13b),  we  find  an  error  of  lower  order  of  magni' 
tude  than  above  when  we  try  to  satisfy  the  body  boundary  condition  above  y  =  0 
This  error  is  not  negligible  even  if  d3R/d63  =0  at  0=0,  and  so  we  assume 
that  this  derivative  does  vanish  there.  However,  we  assume  that  the  curvature 
at  0=0  does  not  vanish,  that  is,  C(0)  f  0  . 


We  again  write  (^(Xjy)  in  terms  of  a  function  of  a  complex  variable,  as 
in  (15)  and  (16)  « 

<j>l  (x,y)  =  Re{  fj(z)}  =  Re{g1(z)}  +i?e{h1(z)}  . 

It  is  also  convenient  to  introduce  the  notation: 

Y \  (x,y)  =  Re{ gj(z)}  ,  (46) 

(x#y)  =  Re{ hjtz)}  .  (47) 

Instead  of  introducing  an  expression  like  (17)  for  hj(z)  ,  we  now  set 

h:  (z)  =  -  j  fj(z){l-  [fo(z)]2}  .  (48) 


Then  we  observe  that 


hl(2)  =  ^lx  " 


j  f”(z){l  -  3  [f  ^  (z)  ]  2 }  . 


Breaking  the  last  expression  into  real  and  imaginary  parts  and  using  properties 
of  f0(z)  already  established,  we  find  that 

/fefhjtx)}  =  iJj1x(x,0)  =  0  ;  (50a) 

-  rm{hj(x)}  =  <|i j  (x,0)  =  J  <*>oXx{l  _3<)>0x}|  =  *lv(x'0)  •  (50b) 

y  y=0  Y 

The  last  result  shows  why  we  introduced  the  expression  in  (48) :  It  takes  care 
of  the  nonhomogeneous  part  of  the  free-surface  condition,  just  as  the  expression 
in  (17)  did.  The  new  form  of  hj(z)  also  has  the  useful  property  given  in  (50a) 
(If  the  form  in  (17)  has  this  property,  it  is  certainly  not  obvious.)  The  other 
part  of  the  complex  potential  satisfies  a  homogeneous  free-surface  condition: 

-  Imlgitx)}  =  y jy (x , 0)  =  0  .  (51) 

The  function  Yj(x*y)  must  be  determined  so  that  (51)  is  satisfied  and  so  that 
the  sum  Yj(x»y)  +li'i(x*y)  satisfies  the  body  boundary  condition  on  <t>j(x,y)  . 

The  contrast  between  (50a)  and  (51)  is  crucial  to  what  follows:  i|>j(x,y) 
is  an  odd  function  of  y  (except  possibly  for  an  additive  constant) ,  whereas 
y^Xjy)  is  even  with  respect  to  y  .  It  is  this  fact  that  leads  to  the  major 
error  in  the  satisfying  of  the  body  boundary  condition,  as  we  now  demonstrate. 

We  should  determine  Yi<x»y)  so  that  it  satisfies  (51)  and  also  so  that 


16 


-r —  =  -  — r —  for  -IT  <  b  <  0  .  (52) 

3n  3n 

Just  as  with  4>g  ,  we  cannot  prescribe  y ^  further,  since  the  determination  of 
y  in  y  <  0  completely  determines  y1  in  y  >  0  .  From  (50a)  and  (51),  it 
then  follows  directly  that 

9y  i  dit>i 

— —  =  +  r —  for  0  <  9  <  IT  (53) 

3n  dn 

So  y1  does  not  cancel  the  effect  of  if>,  in  this  region.  In  fact,  it  pre¬ 
cisely  doubles  it.  Thus  we  find  that 

3(f).  by,  3  4*, 

-r-±-  =  2  =  2  for  0  <  9  .  (54) 

3n  3n  dn 

Since  ip i (x , y )  has  been  completely  specified  already,  we  can  conveniently  use 
it  for  the  calculation  of  3<j>j/3n  on  the  body  above  the  line  y  =  0  .  This 
normal  velocity  component  on  the  body  surface  will  ultimately  have  to  be  can¬ 
celed  by  the  wave-motion  solution. 

The  expression  for  the  normal  velocity  component  on  r  =  R(0)  ,  given  pre¬ 
viously  in  the  form  in  (38),  can  be  rewritten: 

3M  _  „-l-ie  m>  -in’  (8)1 


-sj-  - 


1  r=R(0) 


/rZ  +  R' ' 


We  need  to  evaluate  the  right-hand  side  just  for  small  0  (or  y  ) .  A  few  sub¬ 
sidiary  results  are  useful: 

hj(z)  »  -  ~  {*oxx(xO'°>  +  iy^Oxxx^O'0)  }  '  with  x0  =  R(0)  ;  (56) 

e10  «  1  +  i0  «  1  +  ±2-  ;  (57) 


R(6)  -  iR’ (0)  =  xQ  -  i0R" (0)  =  x0  -  - 

[r2+r.2]-1/2  a  I  . 

*0 

Substituting  all  of  these  back  into  (55)  ,  we  find  that 


iyR" (0) 


1  r=R(0) 


“  I  y^0Xxx(x0,O)  +  c^°^0xx  ' 


One  further  simplification  can  be  made.  From  the  boundary  condition  on 
4>  (see  the  Appendix)  ,  we  have 


Using  this  in  (60)  and  then  finally  returning  to  (54) ,  we  have  our  desired 
result: 


=  2  y  C (0)  <f>0  (x0,O)  for  0  <  0  .  (62) 

r=R(0)  XX 

Since  the  upper  limit  of  the  wetted  body  surface  is  at  y  =  e/2  ,  we  see  that 

the  quantity  in  (62)  is  0(e)  .  However,  it  must  also  be  multiplied  by  e  in 

the  asymptotic  expansion  of  the  full  solution,  and  so  the  error  in  satisfying 

the  body  boundary  condition  is  0(e2)  . 

The  potential  ^(Xjy)  ,  like  <f>Q(x,y)  ,  is  not  analytic  at  (xQ,0)  .  It 
is  interesting  to  compare  their  behavior  in  a  neighborhood  of  this  point.  If 
we  set  £  =  (x-Xg)  + iy  ,  we  find  from  the  Appendix  that  f Q (z)  has  a  term  that 
is  proportional  to  S4log  (,  .  On  the  other  hand,  fj(z)  will  contain  whatever 

singularities  are  contained  in  both  g^z)  and  hj(z)  .  Using  the  result  of 
the  Appendix  in  (48),  we  find  that  the  leading  singular  term  in  hj(z)  is  pro¬ 
portional  to  S3log  C  ,  which  is  still  a  very  weak  singularity.  However, 

g^z)  must  have  a  term  of  the  form  i;2log  t,  .  To  show  this,  we  note  that  (60) 

is  valid  for  y  positive  or  negative.  This  same  expression  gives  Syj/Sn  for 

9  >  0  but  it  gives  the  negative  of  dy-^/dn  for  0  <  0  .  Therefore  3Yj/3n 
must  be  proportional  to  |y|  .  Such  a  property  requires  that  gj(z)  have  the 
form  stated. 


3£l 

3n 


CONCLUSIONS 


There  are  two  principal  conclusions: 

(1)  The  solution  of  the  complete  problem  should  De  broken  into  two 
terms,  as  in  (30)  and  (31),  with  the  first  representing  the  so-called  naive 
expansion.  This  procedure  leads  to  a  homogeneous  free-surface  condition  on 
tile  wave  part  of  the  solution. 

(2)  The  naive  expansion  fails  to  satisfy  the  body  boundary  condition 

on  trie  small  segment  between  y  =  0  ,  the  level  of  the  undisturbed  free  surface, 
and  the  actual  location  of  the  stagnation  point,  which  is  locat.ed  at  y  =  l,  2  . 
Tiie  normal  component  of  velocity  on  the  body,  resulting  from  this  fact,  is 
given  approximately  by  (62) .  The  wave  part  of  the  solution  has  to  cancel 
this  velocity  component. 

The  wave  solution  is  derived  by  Chen  and  Ogilvie  (1982) .  It  may  be  noted 
that  the  form  of  their  solution  differs  in  one  respect  from  that  supposed  in 
(62) :  The  power  of  e  in  the  second  term  of  the  solution  is  initially  taken 

as  unknown,  a  quantity  to  be  determined  during  the  solution  process.  In  fact, 
the  wave  term  turns  out  to  be  higher  order  than  stated  in  (62) .  Its  order  of 
magnitude  is  completely  specified  eventually  by  the  requirement  that  there  be 
a  stagnation  point  on  the  body  at  y  =  e/2  . 
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SINGULARITIES  IN  THE  NAIVE  EXPANSION 


APPENDIX . 


In  the  first  term  of  the  naive  expansion,  4>g(x,y)  ,  there  is  a  stagnation 

point  at  x  =  Xg  =  R(0)  .  This  is  to  be  expected,  of  coarse,  in  view  of  the 
interpretation  of  <J>q  in  terms  of  the  double-body  (rigid-wall)  flow.  it  is 
indicated  explicitly  by  the  combination  of  the  conditions  (11)  and  (13a) .  A 
stagnation-point  flow  is  so  elementary  and  so  well-behaved  that  one  can  easily 
overlook  the  fact  that  it  may  well  be  singular  (in  a  mathematical  sense) .  We 
show  this  explicitly  here. 


We  rewrite  the  body  boundary  condition  on  $  ,  (13a) ,  in  a  form  similar 

to  (55) 


3<Pp 

3n 


r=R 


0 


(0) 


Re 


i0 


fn(z) 


Rq ( 0)  -  iRn(O) 


0 ' 


/r£  +  r’*  2 


0  . 


(Al) 


Since  we  are  considering  the  flow  around  a  double  body,  we  have  taken  the  equa¬ 
tion  describing  the  body  as  in  (34).  Our  purpose  is  to  show  that,  in  order  to 
satisfy  this  condition,  fQ(z)  cannot  be  analytic  at  (xQ,0)  .  Furthermore, 
we  wish  to  determine  the  nature  of  the  nonanalytic  behavior  there. 


Our  approach  is  to  expand  the  quantities  in  (Al)  with  respect  to  the  singu 
lar  point.  The  expansion  of  Rg(9)  has  already  been  given  in  (36),  from  which 
the  expansion  of  Rg(0)  can  be  found  trivially.  Of  course,  there  is  no  diffi¬ 
culty  in  expanding  the  factor  exp(i8)  for  small  0  .  It  is  then  convenient 
to  change  from  0  to  y  by  using  the  relationship 


y  =  Rg ( 0 )  sin  0 


=  xQ0  + 


Putting  all  of  this  together,  we  obtain 


i0  Rn (6)  ~ iR0 (0) 


1  2 

—  ir 


=  1  +  iCy  +  — 


+  •  •  • 

Rg(0) 

are  to 

be 

evaluated 

,2 

f2R'o 

p  ii  2 

_  V . 

] 

iK^ 

+ 

X3 

xo 

X4 

xo 

’  1? 
xo 

1  T3  , 

xo 

(A2) 


+  ...  .  (A3) 


We  note  that  corresponding  to  (A2)  we  can  also  write: 


X  -  Xq  =  Rg  (0)  COS  6  -  X0 

=  -  J  02X2C  +  ~  K03  +  ...  . 

In  these  expressions,  C  is  the  body  curvature  at  0=0 


C  =  C(0)  = 


x0  -  R o  ( 0 ) 


( Cf .  (44).)  We  assume  that  curvature  is  continuous  for  the  double  body. 

We  can  satisfy  the  boundary  condition  if  we  set,  with  5  =  z-x0  , 

fQ(z)  =  A2C2  +  A3S3  +  A4C4  +  ...  +  (log  C  )  [B^4  +  ...]  .  (A6) 

We  write  the  derivative  of  this  expression, 

f  q  (z)  =  2A2£  +  3A  342  +  (4A4+B4K3  +  ...  +  4B4C3logC  +  ...  .  (A7) 

We  evaluate  this  on  the  body,  using  (A2)  and  (A4) ,  and  substitute  into  the  body 

boundary  condition.  We  find  that  the  coefficient  of  each  power  of  y  equals 
zero  if 

1  2K 

A2  =  2  *0Xx(X0'O)  ;  A3  =  "  °  A2  *  B4  =  Tirx3  A2  *  (A8) 

That  is. 


(z)  =  1  'i’Oxx^O'05  {^2  -  c?3  +  A4C3  +  ...  +  C4  log  r,  +  ...}  . 


The  fourth  derivative  of  fQ(z)  is  undefine<^  (infinite)  at  (xQ,0)  .  Thus  we 
have  found  the  leading  singular  term  in  the  double-body-flow  potential. 

From  (A9) ,  we  find  immediately  that 

f;<vi01  ■  Wv°> 


(a  trivial  result) ,  and  also  that 


f  *6*  (x0+i0)  =  -  3  C (0)  *0xx(x0,O) 


Since  the  last  expression  also  equals  ^0Xxx + ^^Oxxy  '  ^Oxxy^X'^  =  ®  from 

(13a) ,  we  have  the  result  needed  in  (61) . 
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